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ABSTRACT 


In  this  report,  an  atteii5)t  is  made  to  determine  svdtable 
smoothing  intervals  to  he  \ised  in  the  smoothing  of  radar  trajectory- 
data.  This  is  done  by  \ising  different  fitting  in-tervals  in  the 
smoothing  of  simulated  noisy  trajectory  data  of  different  signal- to- 
noise  ratios.  The  errors,  or  differences  be-tween  the  true  noiseless 
data  and  the  smoothed  noisy  data,  were  observed  and  plot-ted.  Minimum 
errors  vere  obtained  vhen  the  range  smoothing  interveLL  -was  approximately 
1  to  2  seconds  and  the  elevation  and  azimuth  smoothing  integrals  10 
seconds.  The  plots  suggest  that  elevation  and  azimuth  errors  can  be 
further  reduced  by  smoothing  over  Intervals  longer  than  10  seconds 
but  this  would  not  be  practical  in  many  actvial  missile  shots. 
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INTRODUCTION 


An  attempt  has  heen  made  to  ascertain  appropriate  fitting  intervals 
for  the  smoothing  of  AMRAD-type  trajectory  data.  It  is  felt  that  the 
resxilts  obtained  will  be  \iseful  in  processing  data  from  ATHENA  shots, 
DESCRIPTION 

A  noiseless  trajectory  was  generated  (20  points/sec)  to  simulate 
a  typical  AMRAD  shot.  Different  noise  levels*were  added  to  this  basic 
trajectory  so  as  to  simulate  the  following  signal-to-noise  ratios:  0  db, 

5  db,  10  db,  00  db  (noiseless).  The  resultant  noisy  trajectory  data 
was  then  smoothed,  separately  in  range,  elevation,  and  azimuth,  by 
standard  quadratic  least  sqToares  proced\u*es  of  the  "sliding  arc"  type 
and  the  following  qviantities  were  bbtained:  range,  range  rate,  range 
acceleration  (rate  of  range  rate),  elevation,  elevation  rate,  azimuth, 
azimuth  rate,  and  drag  velocity.  The  corresponding  noiseless  trajectory 
values  were  subtracted  from  these  and  the  differences,  or  errors,  were 
plotted  versus  time.  The  above  procedure  was  repeated  using  different 
smoothing  Intervals  and  the  stated  signal- to-noise  ratios.  Also,  the 
rms  errors  (l.e,,  the  standard  deviation  of  the  differences)  were  computed, 
tabulated,  and,  for  the  5  db  case,  plotted  versus  smoothing  interval. 

Other  parameters  such  as  drag  acceleration,  wei^t-to-drag  ratio,  height, 
path  angle,  etc.,  were  computed  and  plotted  versus  time  for  selected 
cases. 

RESULTS 

In  Figure  1,  the  two  plots  on  the  left  describe  the  geometry  of  the 
♦See  the  Appendix  for  a  description  of  the  noise  generation. 
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simulated  trajectory  giving  latitude  versus  longitude  and  height  versus 
surface  range  with  AMRAD  located  at  the  origin.  No  tangent-plane  approx¬ 
imations  are  involved  in  the  height  vers\as  surface  range  plot.  That  is, 
height  and  surface  range  are  measured  above  and  along  the  sxirface  of  a 
spherical  earth.  The  two  figures  on  the  right  show  range,  elevation, 
and  azimuth  for  the  noiseless  and  S/n  =  5  cCb  cases. 

Figures  2  through  8  show  the  various  errors  (i.e.,  computed  value 
minus  true  value)  in  range,  range  rate,  range  acceleration,  elevation, 
elevation  rate,  azimuth,  and  azimuth  rate  for  the  causes  tested.  Figure  9 
gives  the  drag  velocity  error  versus  time  for  the  cases  tested.  The 
"drag  velocity"  Vj)  is  the  speed  of  the  missile  relative  to  the  radar, 
or  equivalently,  it  is  the  airspeed,  assuming  a  rigid  rotating  atmosphere 
with  no  local  winds. 

Plots  of  the  standard  deviation  (or  rms  values)  of  the  above  errors 
versus  smoothing  interval  for  the  case  of  5  db  signal- to-noise  ratio 
are  given  in  Figures  10  and  11.  Tables  I  throtigh  HI  list  the  rms  errors 
for  all  cases  tested. 

Figure  12  shows  weight-to-drag  ratio,  drsig  acceleration,  drag  velocity, 
height,  path  angle,  and  viewing  (or  aspect)  angle  versus  time  for  the 
5  db  signal- to- noise  ratio  case  and  different  smoothing.  The  path  angle 
GAM  is  the  angle  between  the  velocity  vector  and  the  local  hortzontal. 

The  viewing  (or  aspect)  angle  PSI  is  the  angle  between  velocity  vector 
and  the  line  of  sight.  TFITR,  TFITE,  and  TFITA  are  the  smoothing  Intervals 
for  range,  elevation  and  azimuth  respectively.  The  expressions  vised  for 
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the  computation  of  all  parameters  in  Pig.  12  aie  as  reported  in  47G-6, 
"Determination  of  Weight-to-Drag  Ratio  from  Radar  Meaisxirements,"  8  March  1963^ 
S.  P.  Catalano,  H,  Schneider.  Eqmtion  (33)  of  the  above  report  was  xised 
for  con^mtations  of  weight-to-drag  ratio.  Note  that  weight- to-drag  ratio 
is  best  estimated  in  regions  of  hi^  acceleration. 

CONCLUSIONS 

Referring  to  Tables  I  through  IH  and  considering  the  noiseless 
case,  smallest  rms  errors  were  obtained  when  the  shortest  smoothing 
intervals  were  used.  This  is  as  expected.  With  no  noise,  smoothing  is 
unnecessaiy  as  it  only  corrupts  the  signal.  As  noise  is  eidded,  the 
smoothing  shovild  suppress  the  noise  more  than  the  signal,  thereby  improving 
the  signal-to-noise  ratio.  (This  will  be  true  if  the  noise  fluctuates 
much  more  rapidly  than  the  signal  varies.)  With  excessively  long  smoothing 
intervals,  the  signal  may  also  be  seriously  distorted,  thus  losing  any 
advantage  gained  by  noise  suppression. 

Por  a  signal-to-noise  ratio  of  5  <ib,  minimum  rms  errors  were  obtained 
when  the  range  smoothing  interval  was  approximately  1  to  2  seconds,  and 
the  elevation  and  azimuth  smoothing  intervals  10  seconds.  Pigure  11 
suggests  that  elevation  and  azimuth  errors  can  be  further  reduced  by 
smoothing  over  intervals  longer  than  10  seconds,  but  this  would  no  be 
practical  in  many  actual  missile  shots.  Prcm  Pigure  10  it  is  seen  that 
the  minimum  errors  for  range,  range  rate,  and  range  acceleration  do  not 
occur  at  the  same  range -smoothing  interval.  The  minima,  however,  are 
reasonably  close  together  and  quite  broad,  so  that  very  good  res\ilts 
should  be  obtained  using  a  range  smoothing  interval  from  1  to  2  seconds. 
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A  range  smoothing  interval  vithin  these  limits  should  be  quite  s\iitable 
with  other  signeil-to-noise  ratios  as  may  be  seen  from  Table  I. 

Drag  velocity  and  weight-to-drag  ratio  are  used  as  initial  conditions 
in  a  trajectory  prediction  program.  Figures  9  and  12  should  be  useful 
in  assuring  sufficient  acctiracy  of  these  quantities. 

The  results  in  this  report  should  be  considered  as  a  useful  guide 
rather  than  absolute  rules. 
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APPENDIX 


Discussion  of  Simulated  Trajectoiy  Data  Generation 
The  initial  conditions  at  the  start  of  the  generated  trajectory 
were  as  follows: 

Time  400,000  seconds 


Range 

Range  Rate 
Elevation 
Elevation  Rate 
Azimuth 
Azimuth  Rate 
Weight-to-Drag  Ratio 


499.791  kft 
-20.082  kft/sec 
31,464  degrees 

,016  degrees/second 
-25.00  degrees 
0,0  degrees/second 
110  Ibs/ft^ 


Prom  the  above  conditions,  the  starting  altitude  was  265,160  kft,  emd 
the  drag  velocity  20,082  kft/sec,  A  trajectory  generation  program, 
using  a  spherical  rotating  earth  with  gravity  and  a  realistic  atmosphere, 
calculated  range,  azimuth,  and  elevation  for  the  next  999  points  spaned 
.05  seconds  apart.  At  the  last  calculated  data  x>olnt,  the  range  and 
altitude  were  87.359  kft  and  33.113  kft  respectively.  Gaussian  noise, 
consistent  with  the  desired  slgneQ.-to-nolse  ratio  was  then  added  to 
the  trajectory  data  and  the  resvilting  noisy  data  were  written  on  magnetic 
tape. 


The  standard  deviation  for  the  Gausslem  noise  added  to  the  signeiL 
was  computed  as  described  below.  (These  formulas  are  based  on  the 
assumption  of  a  large  slgnal-to-noise  ratio.  While  this  was  not  always 
the  case  here,  it  was  felt  that  useful  results  would  still  be  obtained). 
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For  range  data* 


'’r-CDsTr  -Ci) 


jLe _ 

2(s/n) 


where  Op  =  ruis  range  error 

Q 

C  =  velocity  of  light  =  9»8351^  x  10  ft/ sec 
6T  =  rms  error  in  estimating  leading  edge  of  p^^lse 
tj.  =  rise  time  of  pulse 
S/n=  signal-to-noise  (power)  ratio 

The  pulse  rise  time  was  assumed  to  he  35  nsec  corresponding  to  a  signal 
bandwidth  of  approximately  20  Mcps.  For  elevation  and  azimuth  data 


=  0.628  9b 
2(s/n) 

where  a  =  rms  elevation  (or  azimuth)  error 

H 

Qg  =  antenna  beamwidth  in  elevation  (or  azimuth) 
s/n  =  signal-to-noise  (power)  ratio 


The  mits  of  cTq  will  be  identical  to  those  of  93.  For  the  simulations 
93  was  taken  as  1*  in  both  azimuth  and  elevation. 

Since  the  noise  bandwidth  of  the  AMRAD  antenna  servo  system  is 
approximately  1  cps,  the  noise  samples  to  be  added  to  azimuth  (or  elevation) 
data  would  not  be  independent  from  point  to  point  with  a  rep  rate  of  20 
point s/sec.  The  following  formula  was  used  to  generate  the  azimuth  (and 
elevation)  noise. 

♦Skolnik,  M.  I.,  "introduction  to  Radar  Systems^"  McGraw-Hill  Book  Company, 

1962,  pp.  462-464. 

**Skolnik,  op.  cit.,  pp.  4t6-4T7. 
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j»k+N-l 

■v '  I"j 

J  =k 


th 


=  Gaussian  noise  to  be  added  to  k  azimuth  (or  elevation) 
sample 


,th 


m.  c  J  independent  Gaussian  noise  san5)le  of  standard  deviation  cr* 

J  9 


/Repetition  Rate  \  l/2 
”  ^Koise  Bandwidth-' 


This  assvired  that  the  noise  would  be  correlated  over  a  1  second  interval. 
The  calculated  and  observed  rms  noise  values  are  given  below: 
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Slgnsil-to-Nolse  Ratio 

Range  Smoothing  Interval 
(seconds) 

RMS 

Range  Error 
(feet) 

RMS 

Range-Rate  Error 
(ft/ sec) 

RMS 

Range  Acceleration  Error 
(ft/sec2) 

Noiseless 

.5 

2.0 

1.1 

3.6 

2.0 

1.9 

ll^.I^ 

5.2 

^^.0 

3.7 

55.3 

15.5 

10  dh 

2.0 

2.1 

ll^.5 

6.8 

5  db 

.5 

h,2 

13.8 

194.6 

1.0 

3.^ 

6.4 

39.4 

1.5 

3.1 

8.8 

l4.i 

2.0 

3.0 

l4.6 

8.6 

h.O 

^^.3 

55.4 

15.6 

0  db 

2.0 

3.3 

15.9 

13.4 

TABIE  I 

RMS  Errors  of  Range,  Range  Rate,  and  Range 
Acceleration  for  Different  Range  Smoothing 
Intervals  and  Slgnal-to-Nolse  Ratios 


Signal-to-Noise  Ratio 

Elevation  Smoothing  Interval 
(seconds) 

RMS 

Elevation  Error 
(degrees) 

R>© 

Elevation  Rate  Error 
(degrees) 

Noiseless 

2.0 

.001 

.0003 

5.0 

.0002 

.002 

10.0 

.003 

.008 

10  db 

2.0 

.097 

.092 

10.0 

.045 

.012 

5  db 

2.0 

.184 

.176 

5.0 

.135 

.056 

10.0 

.108 

.028 

0  db 

2.0 

.413 

.384 

10.0 

.289 

.051 

TABIE  n 

RMS  Errors  of  Elevation  and  Elevation  Rate 
for  Different  Elevation  Smoothing  Intejrvals 
and  Slgnal-to-Nolse  Ratios 


H 

O 


Slgnal-to-Noise  Ratio 

Azimuth  Staoothing  Interval 
( seconds ) 

RMS 

Azimuth  Error 
(degrees) 

RMS 

Azimuth  Rate  Error 
(degrees) 

Noiseless 

2.0 

.00001 

.00001 

5.0 

.00001 

.00001 

10.0 

.00002 

.00004 

10  dh 

2.0 

.102 

.110 

10.0 

.o4o 

.010 

5  dh 

2.0 

.194 

.198 

5.0 

.138 

.069 

10.0 

.091 

.019 

0  dh 

2.0 

.423 

.453 

10.0 

.187 

.o44 

TABIE  m 

RMS  Errors  of  Azimuth  and  Azimuth  Rate 
for  Different  Elevation  Smoothing  Intervals 
and  Slgn€d.-to-Nol8e  Ratios 
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Pig.  4,  Range  acceleration  error  vs,  time  for  different  sigaal-to-noise 
ratios  and  smoothing  intervals. 
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Pig,  5.  Elevation  error  vs.  time  for  different  signal- to-nolse  ratios 
and  smoothing  Intervals. 
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Fig.  5  continued 
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Pig.  6,  Elevation  rate  error  vs,  time  for  different  signal- to-nolse 
ratios  and  smoothing  Intervals. 
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Fig.  6  continued 
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Pig,  T.  Azinnith  error  vs,  time  for  different  signal- to-nolse  ratios 
and  smoothing  Intervals, 
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Fig.  7  continued 
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Fig,  8,  Azlaruth  rate  error  vs,  tine  for  different  signal- to-nolse 
ratios  and  smoothing  Intezvals, 
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Fig.  8  continued 


Pig,  9,  Drag  velocity  error  vs,  time  for  different  smoothing  intervals 
and  signal- to-noise  ratios, 
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Fig.  9  continued 
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Pig.  11.  RMS  angle  errors  vs,  smoothing  interval  for  the  S/N  ■  5  case. 
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Pig.  12.  Trajectory  parameters  vs.  time  for  the  S/n  »  5  db  case  with 
dlffeirent  smoothing. 
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